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Abstract We obtain an approximate solution for the motion of a charged particle 
around a Schwarzschild black hole immersed in a weak dipolar magnetic field. We 
focus on eccentric bound orbits in the equatorial plane of the Schwarzschild black 
hole and derive an analytic expression for the spectral distribution of the electro¬ 
magnetic emission from a charged particle on such an orbit. Two sets of harmonic 
contributions appear, with specific frequency spacing. The expression can be written 
in compact form, if it is truncated up to the lowest order harmonic contributions. 


1 Introduction 

The origin of astrophysical magnetic fields on galactic and extragalactic scales re¬ 
mains one of the open questions of modern astrophysics. Microgauss magnetic fields 
have been observed in elliptical and spiral galaxies, in the intracluster medium, and 
in damped Lyman alpha systems over cosmological distances (see e.g. [1] for a re¬ 
view). In the past, several investigations of bremsstrahlung radiation were carried out 
in the case of a Schwarzschild black hole embedded in a large-scale homogeneous 
asymptotically uniform external magnetic field ([2-4]). The total radiated power was 
obtained through a generalized Larmor formula based on arguments about covari¬ 
ance under Lorentz transformations ([5, 6]). In a broad sense, bremsstrahlung [7] is 
the radiation emitted by an accelerated charged particle. The angular distribution of 
the radiation produced from an orbiting charged particle in the presence of a uniform 
magnetic field has been discussed extensively with the aid of the Lienard-Wiechert 
potentials at large distances in [6, 8, 9]. Along the way, several interesting features of 
the dynamics of charged particles in the presence of both the black hole gravitational 
field and the uniform magnetic field were investigated. In [10] and references therein, 
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the authors explored two different types of orbits, with and without curls. This anal¬ 
ysis was generalized by [11] in the case of a Schwarzschild black hole embedded 
in a dipole magnetic field. An interesting result obtained in [12] was the possibil¬ 
ity of trapping relativistic charged particles in a potential well inside the position of 
the innermost stable circular orbit. In the present work, we investigate the eccentric 
orbits of charged particles around a Schwarzschild black hole embedded in a dipole 
field and obtain their electromagnetic radiation spectrum, assuming slow motions and 
weak fields. Our paper is organized as follows; In Section 2 we present the mathemat¬ 
ical formalism. In Section 3, we investigate eccentric bound orbits in the equatorial 
plane and in Section 4, we discuss the angular distribution of the electromagnetic 
radiation from those particles. We close with a discussion in Section 5. 


2 Equations of motion in Schwarzschild space-time 


We follow the 3-1-1 formalism of [5] and [13] for a dipole magnetic field in a 
Schwarzschild background and work in geometric units with c = G = 1 (in some 
equations the speed of light c is kept on purpose). Greek spacetime indices take val¬ 
ues from 0 to 3 and Latin indices from 1 to 3. In a (f,r, 0,^) coordinate system the 
Schwarzschild metric is 

ds^ = — ^1 — —^ dt^ + rs [dO^ + sin^ 9d^^), ( 1 ) 


where rs = 2M is the Schwarzschild radius of the black hole. The components of the 
covariant four-velocity are = {iP,u'), with = —1. The equations of motion 
of a charged test particle in the presence of a dipole magnetic field can be derived 
from the Lagrangian 


2Jf = 




dx^ q dx^ 

j-j- '' -^ 

dT dx m dx 




( 2 ) 


where x is the proper time and m and q are the particle’s mass and electric charge 
respectively. The Lagrangian now reads 


2j§f = - 




+ r^{e 


' sin^00^) -b — -j—Au. 

^ m dx ^ 


(3) 


where a dot denotes differentiation with respect to x. 

The motion of a charged particle around a Schwarzschild black hole immersed 
in a dipole magnetic field [11], [12] is characterized by two conserved quantities, the 
energy E >0 and the generalized angular momentum L, given by the expressions 


'=£(■- 7 )^ ■ 

. L q As 

^ ^ 2 o 7 ^ 2 o ’ 

sin 0 m sm 0 


(4) 

(5) 


(6) 
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We assume a dipole electromagnetic four-potential of the form = (0,0,0,A,|,), 
where 

A^i, = ew{r), (7) 

and 

W(r) = ln(l-^)+^ + ^. (8) 

Above, jj.ti is the magnetic dipole moment for an observer at infinity. Thus, for the 
metric (1) the normalization of the 4-velocity yields 




= E^ 



l+r^e^ + r^sin^O 




(9) 


We make the following changes in the above equation. First, we change the T-derivatives 
to 0-derivatives, through the relation r = {dr/d(j))(j)', second, we assume that the par¬ 
ticle motion is confined in the equatorial plane (0 = 7t/2); finally, we assume that 
r ^ rg, thus W{r) fv —r|/3r^, and we arrive at the following equation of motion, 
which describes the geometry of charged particle orbits in the invariant equatorial 
plane 


dll 

d^ 


E^-l 


+ ^u-u^+{rs-2El)i. 


{l+Qu)^[ L2 ' i2 

+ {2Qrs-n^)u‘^ + n^rsu^ 


( 10 ) 


where ii^Xjr and 



( 11 ) 


The parameter Q. can be negative or positive, depending on the sign of the charge 
q and on the azimuthal generalized angular momentum L. In what follows, we will 
restrict ourselves to positive values of Q.. 


3 Bound orbits 

The motion of charged particles in a dipolar magnetic field around a static black hole 
was studied in [11, 12] and [15] (see also the review in [2]). A synchrotron mechanism 
was considered for circular orbits of charged particles in the presence of a uniform 
magnetic field. Here, we focus only on eccentric bound orbits of charged particles of 
the form r = l/{\ -becos^), where e is the eccentricity of the orbit with 0 < e < 1, 
I is its latus rectum and ^ € [0,27r] in the equatorial plane of a Schwarzschild black 
hole, immersed in a dipole magnetic field, when r/rg ^ 1. We seek approximate, 
analytic solutions, by keeping only orders up to 0{ir') in the series appearing in the 
nominator of (10). With theses assumptions, we find values of £, L and Q such that 
Eq. (10) admits bound orbits. More details on this method can be found in [16-19]. 
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The geometry of the orbits is determined by the real roots of the equation f\(u) = 0 
where 

/ i ( m ) =- -^ + -^ii-i? + {rs-2a)u\ (12) 

(when l+£2u^0). The orbital eccentricity e and latus rectum I = a(l — e^) (with a 
being the major axis of the ellipse) satisfy 


g = i[2Z-(3 + .2)(r,-2f2)], 
il^ = l[/-2(r,-2f2)](l-.2). 

In this case, the equation /i («) = 0 admits three real roots which are 


1 —e 


Ml = 


1 +e 


1 


-Y~ ’ “2 = —^ , M3 = -^ - 7 

I I rs — 2L2 I 


From (13) we deduce that 


M 1 Q. 




Z - 3+e2 ' ; ’ 

and assuming that M 3 > M 2 we arrive at another interesting inequality 

I 


rg — 20. < 


3 + e 


Substituting (13), (14) and 

into ( 10 ), we obtain 
'd^'' ^ 


: = y(l +ecos<*), 


1 


J [1 + (f 2 /Z)(l+ecos^)]2 

which reduces to 


1 — (3 +ecos^) 


rs-20 

I 


O 

l + -(l-e) 


d^ 


20e\ cos'^ j2)d^ 


= ±\fXd(^, 

with A = 1 — (3 — e) > 0 and 

2e{rs-20) 


k^ = 


l-i3-e){rs-20) 


< 1 . 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 


Integrating (20) we find 


2[Z + f2(l-e)]F[cos(<s/2),;fc]-4f2e{F[cos(^/2),Z:]-£[cos(^/2),Z:]} = ±VAZ0, 

( 22 ) 
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where F[cos (i*/2),A:] and ^[cos {^ /2), k] are the elliptic integrals of first and second 
kind, respectively, with modulus k <1. 

Expanding the elliptic integrals in terms of the small modulus k and keeping only 
zero order terms in k < 1 (see [20-22]) , we obtain 

Qe 


Itoii-I.) < 1 mil j; = rdfraj- 
From (23) we compute cos^, which appears in (18) and obtain 

cos^ = - cos 0 ) T ^ Sin^ (x 0 ). 

Substitution of (24) into (18) yields 

. = 1 [l - .cos(j:« IF (rt) 

To linear order in e, the last equation becomes 





-ecos(x(p)], 


(24) 


(25) 


(26) 


where the origin of (j) may be determined by Eq.(23) setting ^ = 0. 

The approximate solution may now be completed by direct integration of the 
equations 


d(j) 

dx 

dt 

dx 


= Lu^( \ + £lu) , 
= E{\-rsuY^ . 


From (26) and (27) we obtain 




XL 

= TXljdt, 


[l-ecos(:j;0)]2[pi+p2Cos(x0)] Efi 


where we have neglected the 6>(e^) term 2e^(^) and setpi = (l—2p)(l + y) and 
P2=2p-^(l-4p). 

Integrating (28), we find 

Nt = B\ arctan 

+ -83 


(27) 


(28) 


1 —e 
sin ix^) 


1 +e (y0) \ 

tan-^;— +82 arctan 


Pi+gp2 U0) 
Pi—ep 2 2 


2(1 -ecos {X(l>)) ’ 


( 29 ) 
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where 


and 




Bi 


B 3 


2[(pl+p2)+p2(l-e^)] 

(l-e2)3/2(pi+p2)2 

2p| 

{pi+p2)^\/pi-e^p2 

2e 

il-e^){pi+P2)’ 


_ XB L \/l-(3-e)^ 

1 + ^(1 _ 2 e ) 


(30) 


(31) 


has dimensions of frequency. In Appendix B, we compute the right-hand of (31) in 
terms of e, p and Q ^0 with the aid of (13), (14) and find a rather lengthy expression. 
However, for = 0, we compute the ratio L/E using (13) and recover the Newtonian 
frequency of a Keplerian orbit with the same e and /. 

Note that, when both the eccentricity e £ (0,1) and p =M/l = rs/{2l) are small, 
the coefficients Bi, B 2 and Bj can be approximated as 


Bi=2 + 0(e)^2, B 2 = 0(p^), B^ = 0{e). (32) 


Thus, (29) takes the approximate form 


tan 




tan(Af), 


(33) 


since the second and the third term are very small in comparison to the first one 
multiplied by Bi. Furthermore, using (33) we find 


cos{x(j)) = 


cos {2Nt) — e 


1 — ecos {2Nt) ’ 

,in(rt)=v£^( 

\ 1 — ecos {2Nt)) 


(34) 


Finally, because of our previous approximations in (32), Eq. (26) yields 


/[I — ecos(2Af)] 

(1 +e2) —2ecos{2Nt) ' 


r = 


(35) 
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4 Spectral distribution of radiation 


We are next interested in examining the electromagnetic emission of a particle with 
charge q that is following an eccentric bound orbit, such as the one discussed in 
Section III. Our aim is to obtain analytic expressions for the distribution of energy 
as a function of angle and radiation frequency (see e.g. [6, 8, 9]), valid under the 
assumptions made in Section III. We make use of the well-known result [6, 8, 9] 


d^<g 

dco'dQs 


q^co'^ 


y” n X (n X 


(36) 


where S is the total energy radiated per unit frequency O)' per unit solid angle 12^ [8] 
in which we replace the internal and external vector products with general-relativistic 
expressions (derived in the 3 4-1 formalism [12, 13]). For our further derivations, it 
is necessary to know the spatial velocity v and position r(f) over a small arc of the 
trajectory whose tangent is pointing towards the observation point. In our case the 
segment of the trajectory lies in the .ry-plane, and the unit vector n in the direction of 
radiation lies (without loss of generality) in the yz-plane at an angle X with the z-axis. 
In this case, the unit vector n is 


n = sinACy-fcosACj, (37) 

where e,, and Cj are the unit vectors along the y and z axes, respectively. In [6, 8] the 
authors examined the angular distribution of the energy radiated per unit frequency 
cu' per unit solid angle for X k, nj2, since appreciable radiation intensity is con¬ 
centrated mainly in the equatorial plane of the orbit, and the duration of the pulse is 
short. Here, we find an expression for any angle X G (0,2n) (but which is valid only 
under the assumptions mentioned in Section III). 

The internal and external vector products in (36) are derived in spherical polar 
coordinates and become * 

■ • r 

n • r = Yijn'r^ = — [sin A sin 9 sin0 -f cos 0 cos A], (38) 

^J^-ralr 


and 


n X (n X v) = Xtr 


-ee- 


rsin 0 


(39) 


where the explicit expressions of X,Y and Z are given in Appendix A. In the equato¬ 
rial plane (0 = ?r/2), (36) then becomes; 


d^^ 

dco'dQs 


q^CO'^ 



Fq Zo 

XoCr H -eg H - e,j, 

im't [ ■ / Tsinlsin^ 1 

e'®* X exp< -10)' - — \dt 

J-oa 

r r ^ 



( 40 ) 


* The cross vector product is (n x v)/ = with £123 = 

the spatial part of the metric. 


■ n 

I- Sin 6 
>^{l-rs/r) ’ 


while in the dot product Yij is 
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where 


^0 = - 


\/(l-r5/r) 


2. dr . 2 T ^d{rcos^) 

cos A— — Sin A cos 0-;- 

dt dt 


Yq = —r sin A cos A 

Z{) = 


d{r sin (j)) 
dt 


. . ,d{rcos(j)) 2 i f d<j)\ 

sin A sin 0-cos A r — 

dt \ dt J 


Using (67), the exponential function in (40) becomes 


^ £ j^^Q^yico'-2mN) ^ 


where 7„, are the Bessel functions of the first kind, while we set 


4mN) 


Qi = cu'sinA- 


l 


and 22 = ®^ sin A 


Eq. (40) with the help of (64) and (44) now reads 


l{2-5n)e 

2(l-2jU)3/2' 


dco'dQs 47t:2 


[(lUr)^ + (lUf)2 + (lk3'")2] 


where 


(41) 

(42) 

(43) 


(44) 


(45) 


(46) 


W{^= f = 2/,W ( 52 d{co'-4Nm)W^l 

J V ^ \n}=—oo 

- £ 5{(0'-lNm)wA , (47) 

m=—oo J 

wm = r ^ -.f^lNsinlX I f, S {co’ - 2 mN)WS, 

+ f;5(a)'-4miV)lU2l^ ,(48) 

w"' = ( e'®'fi-n-r(f)) ^ £ 5(0)' - 2mN)WZ 

J —oo r y Z \ 

+ £ 5 ( 0 '- 4mN)Wil ) . (49) 
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Furthermore, (46) with the aid of (47), (47) and (49) yield 






dco'dQs 471^ 

xj £ 5(co'-2Nm)[{Wi:f + {WSf + (WZf] 


+ ^ 8{(0'-4Nm)[{Wll)^ + {wrt,f + {wri,f^ 


(50) 


We have used the approximation that 5^(0)') = 5(a)') in the limit T ^ where 

T is the radiation emission time during the orbital motion of the charged particle 
around the black hole, while we set 

W^a = 2cos 21(1 -2Ai)i/‘^c[7,„+i(ei) -y„,_i(ei)] 

(l-2Ai)2 


— sin^ A 


■ [Jm+2{Q\) —Jm-liQl)] 


+ 


-2(i-2M)^e[y,„+i(ei)-y,„_i(ei)] 

(1-2m)(7-6m)c 


[Jm+i{Q\) — Jm-iiQl) + Jm+\{Ql) — Jm-\{Q\)] ^ ,(51) 


= COs 2A(1 -2M)‘/^c[72m+l(e2) -^2,„-!(&)] 

^ [Jlm+liQl) - Jlm-liQl)] 


— sin^ A 

— 2(1 — 2^)^e[y2m+l(22) — Jlm-xiQl)] 
(1-2m)(7-6m)c 


+ 


8 


[A2m+3(22) -72m-3(G2) +-/2m+l(22) -72m-l(22)] 2 , 


WZ = [Jm-l{Qx)+J,n+l{Ql)]+2e[J^^2{Qx)+Jm+2{Qx)]. 

^2b = [J2m-\{Q2) + J2m+l{Q2)]+'^e[J2m-2{Q2) + J2m+2{Q2)], 


(52) 

(53) 

(54) 


= sin"A |-[27„,(ei) + y,„_ 2 (ei) - A«+ 2 (ei)] 

5e 

[Am— 1 (2l ) 4" Am+1 (2l ) Am—3(2 i ) Jm+2 (2l )] 

+2ecos^A[7m(Gi) +Am-i(2i) +Am+i(2i)]5 (55) 

5e 

+ -^[A2m-l(22) +A2m+l(G2) -A2m-3(22) - A2m+3(22)] 

+ -COS^ A[4ym(22) +A2m-l(G2) +A2m+l(22)A2m+3(22) — A2m-3(22)]- (56) 

Recall, that the argument Q 2 depends on the geometry of the orbit, since it is propor¬ 
tional to the eccentricity. The analytic expression (50) is our final result. 








10 


D. B. Papadopoulos et al. 


5 Discussion 

We obtained an analytic expression for the spectral distribution of the electromag¬ 
netic emission from charged particles in an eccentric, bound equatorial orbit around 
a Schwarzschild black hole, embedded in a dipolar magnetic field. Our final result 
(50) includes an infinite number of harmonics. We note that in the second sum, the 
argument Qi in the Bessel functions depends on the eccentricity of the charged par¬ 
ticle, but in the first sum, the argument Q\ does not. It is evident from the analytic 
expression (50) that the spectrum of the emitted radiation is composed of spectral 
lines appearing at frequencies intervals Acoi = 2N and A ( 02 = 4A^, where N is given 
by (31). 

Our main result (50) can be further integrated over all frequencies o' and all 
angles, yielding a distribution in terms of the harmonics m. The relative contribution 
of individual harmonics will depend on the precise orbit. For small eccentricities, 
such as those considered here, we expect only the first few harmonics to be important. 
Under this assumption, one could thus truncate the infinite series in (50), in order to 
arrive at a more compact analytic expression. We can obtain a rough estimate of the 
total radiated power for an electron in orbit around a ten solar mass Schwarzschild 
black hole with M = 2 x 10^4 g, 

rs = 3 X 10^ cm, immersed in a dipole magnetic field 
typical of intergalactic space of order B ~ 1 /rG. We further restrict our analysis to 
equatorial orbits with eccentricity e = 0.5 and latus rectum on the order of I ^ risco- 
In that case, L = 0.53 x lO^cm^/sec, ^ ~ 0.81 x 10^ cm, ^ = rs/2l = 1 / 6 , and 
^ ~ 1.41 X 10^^ is indeed less than 1, so our approximations are valid. We obtain 
an estimate of the ratio L/El^ ~ ±0.79 x lO^sec/cm^, and Nl/c ~ 0.141. The total 
radiated power 3^ is then obtained by integrating Eq. (50) over all frequencies and 
solid angles, and by dividing with T, namely ~ 1.4 x 10^^^ erg/sec. This yields 
a timescale for orbit evolution on the order of T ^ nieV^ j^ ~ 10 *® sec, which is 
dynamically unimportant. In a future publication, we will relax our current working 
assumption that ^ ^ 1 , and will also consider higher values of the magnetic field for 
which radiation effects will become astrophysically significant. 
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Appendix A 

The explicit expressions for A, y,Z are 



(57) 

(58) 

(59) 


where the components n\ v' are computed in the curved space time and £123 = 

2 • n 

—p=S£2=. Using (34), (35) and their derivatives, we find; 

^JA-rs/r) 
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Xo = 2 Wsin( 2 M){cosUe(l-e^) x [l-2^(1+ 6 ^) - (1-4^)ecos(2M)] 

X [1 — ecos(2M)]^^^^ [(1 +e^) — 2ecos(2M)] 

— sin^ A(1 —e^) [cos (2Nt ) — e\ 

x[l -2^(1 +6^) - (1 -4^)ecos(2M)]‘^^] 

X [1 -ecos(2M)]^^^^[(l + 6 ^) -2ecos(2M)]^^} , (60) 


-Fo = — 2 W sin A cos A 
r 


a/( 1 — e^)[(l +e^)cos ( 2 M) — 2 e] 
[( 1 + 6 ^) — 2 ecos( 2 M )]2 


(61) 



-IIN 



(I_e2)3/2gijj2 

[1 —ecos( 2 M)][(l +e^) — 2 ecos (2Nt)]^ 


+ COS^/l 


Vl 


[(1 +e^) —2ecos {2Nt)\ 


(62) 


while the term jn ( 40 ) becomes 

rsinAsin 0 r -r , , > r , M-3/2 

- = /v 1 — e2sinAsin(2M) x [1 —ecos(2A^f)] ' 

\/(l-rs/r) 

X [l-2^(l+e2)-(l-4^)ecos(2M)]^‘^^ (63) 

Keeping terms up to e^, an approximate form of the (60-63) is 

Xo = 2Wsin(2M)cos2l[(l -2^)'/'‘e + 0(e^)], (64) 

-Fq = —2Wsinlcosl x {cos(2M) -|-2[cos^ (2A^0 ~ l]^ + 61(e^)}, (65) 

r 

-Zq = —2/N{sin^ A sin^ (2Nt) x [1 -|-5cos(2M)e-|-(9(e^)}, ( 66 ) 

r 


and 


1 


(2-5At) 

(1-2ai)3/2 


cos {2Nt)e + 0{e'^) 


n-r = lsinAsin(2A^f) 


(67) 
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Appendix B 

In (31) we computed the ratio L/E with the aid of (13), (14). In this case, we hnd 


L2 rsl 
Furthermore, (68) gives 


, Q 

1 _( 3 +, 2 )(^-_) 


1 

P 




1_4(^--)(1-.^) 




- =±V^Ul-2M)^-4MV + -[(3+e2)-4^(l-e2)] 


- 1/2 


( 68 ) 


(69) 


In (69), setting = 0 and e = 0, we recover Brener’s result [24]. In the case where 
= 0 and e ^ 0, we recover the result obtained in [16]. 

From (31) and (68) we hnd 


2N = ±'- 


pTT 


I (1 - 2^)2 - + y [(3 + .2) - 4 m(1 - e^)] | 



r , ,1 

-1 

£2 

X 

l + _(l_2c) 

X 

l-(3-e)(2M-2y) 


1/2 


(70) 


In the limit = 0, we recover the coefficient of the Eq. (136) in [16]. Observe, that 
on the right side of (70) the coefficient ±/W*/2//3/2 corresponds to the Newtonian 
frequency of a Keplerian orbit atr = Hnc = G= 1 units. 
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